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g_^ . Abstract 

r^ I . In the modern formulation of lattice gauge-fixing, the gauge fixing condition is written 

in terms of the minima or stationary points (collectively called solutions) of a gauge-fixing 

fT^ \ functional. Due to the non-linearity of this functional, it usually has many solutions called 

^ . Gribov copies. The dependence of the number of Gribov copies, n\U\ on the different gauge 

r^ ' orbits plays an important role in constructing the Faddeev-Popov procedure and hence in 

\l , realising the BRST symmetry on the lattice. Here, we initiate a study of counting n\U\ for 

different orbits using three complimentary methods: 1. analytical results in lower dimen- 
sions, and some lower bounds on n\U\ in higher dimensions, 2. the numerical polynomial 
CO , homotopy continuation method, which numerically finds all Gribov copies for a given orbit 

(«_^ ' for small lattices, and 3. numerical minimisation ("brute force"), which finds many dis- 

Cn ' tinct Gribov copies, but not necessarily all. Because n for the coset SU(A'^c)/U(l) of an 

" . SU(Ai'c) theory is orbit-independent, we concentrate on the residual compact U(l) case in 

^ ' this article and establish that n is orbit-dependent for the minimal lattice Landau gauge and 

• '^ , orbit-independent for the absolute lattice Landau gauge. We also observe that contrary to 

|r% . a previous claim, n is not exponentially suppressed for the recently proposed stereographic 

3 ' lattice Landau gauge compared to the naive gauge in more than one dimension. 

1 Introduction 

Though quantum field theory with the perturbative approach has been extremely successful, 
there are many important physical phenomena, such as quark confinement and dynamical chi- 
ral symmetry breaking in quantum chromodynamics (QCD), which need a non-perturbative 
treatment. To understand such a non-perturbative phenomenon, a genuine non-perturbative 
approach to QCD is essential. An immensely useful approach to study such non-perturbative 
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phenomena is lattice field theory [I] [H [3] • In this approach, the Euclidean space-time is discre- 
tised, so the four-dimensional space-time integral is replaced by a discrete sum over all lattice 
points while derivatives are replaced by finite differences. Then, Monte Carlo methods used in 
statistical mechanics can be applied to calculate the expectation values of observables. 

In the continuum, a promising approach to study non-perturbative phenomena in QCD is to 
study truncated systems of Dyson-Schwinger equations (DSEs) |4j, which are the equations of 
motion for QCD Green's functions. Since each gauge configuration comes with infinitely many 
equivalent physical copies, the set of which is called a gauge-orbit of the gauge configuration, the 
generating functional and hence any relevant quantity requires gauge- fixing to remove redundant 
degrees of gauge freedom. 

The DSE approach has a clear advantage in the low momentum region of QCD. However, lattice 
QCD provides an opportunity to do first principles calculations of non-perturbative quantities 
in QCD: the approximations involved in lattice QCD can be systematically removed, unlike 
the truncations of DSEs. Thus, lattice simulations can provide an independent check on the 
results obtained in the DSE approach. A gauge field theory put on the lattice is manifestly 
gauge invariant, i.e., one does not need to fix a gauge on the lattice to calculate gauge invariant 
observables. However, to compare with DSE results, gauge fixing is necessary, and it is mainly 
for this reason that lattice Landau gauge studies have gained a large amount of interest recently. 

The standard approach of gauge-fixing in the perturbative limit is the Faddeev-Popov (FP) 
procedure [5]. There, a gauge-fixing device which is called the gauge- fixing partition function, 
Zgfi is formulated. With an ideal gauge- fixing condition, Zgf is equal to one. This unity 
is inserted in the measure of the generating functional so that the redundant degrees of free- 
dom are removed after appropriate integration. The generalisation of the FP procedure is the 
Becchi-Rouet-Stora-Tyutin (BRST) formulation [51. The assumption that the gauge- fixing 
condition is ideal, i.e., the gauge-fixing condition has a unique solution, is crucial here. How- 
ever, V. N. Gribov, in 1978, found that in non-Abelian gauge theories a generalised Landau 
gauge-fixing condition treated non-perturbatively would have multiple solutions, called Gribov 
or Gribov-Singer copies [ZllHlll], the effects of which should be properly taken into account. 

There is one more obstacle on the lattice. If a procedure analogous to the FP procedure is carried 
out on the lattice, the lattice analogue of Zgf turns out to be zero [HI HH] due to a perfect 
cancellation among Gribov copieo Thus, the expectation value of a gauge-fixed observable 
turns out to be of the indeterminate form 0/0. This problem is known as the Neuberger 0/0 
problem. In other words, BRST formulations on the lattice can not be constructed. This severely 
hampers any comparison with continuum DSE studies [Ml [151 IE] ■ 

On the lattice, Landau gauge fixing is usually formulated as a functional minimisation problem. 
That is, instead of solving the lattice counterpart of the gauge-fixing conditions, one numerically 
minimises a gauge-fixing functional, whose first derivative with respect to gauge transformation 
is the lattice counterpart of the gauge-fixing condition. The space of minima, called the first 
Gribov region, also has many Gribov copies but there is no cancellation among these Gribov 
copies hence the Neuberger 0/0 problem is avoided. However, the number of minima may be 
different for different gauge-orbits making the corresponding Zgf orbit-dependent. In this case, 
inserting Zgf in the partition function becomes quite cumbersome. It is this orbit-dependence 
or orbit-independence of the Zgf in various versions of lattice Landau gauge which is the main 



^See Refs. |lllll2|[T3l for earlier attempts 



focus of this paper. 

Gribov copies not only play an important role in constructing the BRST symmetry on the lattice, 
but it is also extensively argued that they influence the infrared behaviour of the gauge depen- 
dent propagators of gauge theories (See, e.g., [171 [HI [H], and [20l [2TJ [22j [23] for a continuum 
perspective) . 

This paper is organised as follows. In Section [2j we introduce gauge-fixing on the lattice, fix the 
notation and set up the general problem. We also introduce several versions of lattice Landau 
gauge which we will study in the remainder of the paper. In Section [3] we formulate the gauge 
fixing problem for the compact U(l) case. We also summarise the previously known analytical 
results in one dimension and work out some analytical results in higher dimensions. In Section |31 
we introduce a novel numerical method which finds all the stationary points of a given potential 
provided that the corresponding stationary equations are polynomial-like and have only finitely 
many solutions. We show how to translate the lattice Landau gauge-fixing equations for the 
compact C/(l) case into a system of polynomial equations, and then find all the Gribov copies 
for small but non-trivial lattices. Section [S] contains the results from our numerical minimisation 
of the gauge fixing functional. Finally, in Section [51 we summarise our main results and discuss 
the implications of our findings. 

2 Gauge-fixing on the lattice 

A gauge field theory can be studied nonperturbatively by discretising the Euclidean space-time 
and putting the gauge and matter fields of the theory on a four-dimensional space-time grid [1] . 
The gauge fields are usually defined through link variables t/^^^ £ G where the discrete variable 
i denotes the site index, /it is a directional index and G is the corresponding group of the theory. 
The relation of the gauge-fields on the lattice Ui,^ = U^{x) to their continuum counterpart 
Af,{x), is given by 



where V denotes path ordering. 

The lattice expectation value of an observable O is 
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where S\U] is the lattice action. 

The action S\U] is invariant under Ui,f_, — s- U'l — gJUi^^gi+f,, where the gauge transformations 
gi ^ G are defined at each lattice site. In other words, a set of randomly chosen {U} comes 
with infinitely many physically equivalent configurations the set of which is called a gauge orbit 
and represented by the set {[/}. This gauge invariance makes the integration in the continuum 
counterpart of Eq. ^ ill-defined. In the continuum, to get rid of this ambiguity, one can choose 
to take exactly one representative (or more than one but a finite number) of each gauge-orbit. 
This is called gauge fixing. 



Choosing the representatives of each gauge-orbit can in practice be done by imposing a constraint 
of the form f{U^) = in Eq. ([2]). One takes the integral over the gauge orbit, U^, 

n[U] = Jdg detMppig) 5{f{m)) =: Zgf- (3) 

where Zqf is called the gauge fixing partition function, and det Mpp is the determinant of 
the Jacobian of the gauge-fixing condition, called the Faddeev-Popov (FP) determinant. Here, 
n\JJ\ may be a function of the gauge orbit {C/}, that is, n\U\ = n[C/^]. For non-Abelian theories 
treated nonperturbatively, the surface f(U^) = will in general have more than one solution, 
called Grihov copies. If n\U] is non- vanishing for all orbits, then inserting Eq.(|31) in ^ gives, 
after factoring out the gauge volume, 

/ n.,^ dU^^^^det Mpp{g) S{f{US)) exp[^S[U]]0[U.^, 



10) = '^ ^""'"' - - ^__^_!Z. (4) 

provided that the operator 0[U] is gauge invariant. However, this expression can also be used 
for gauge dependent operators, in which case it defines the gauge fixed expectation value. 

The most popular gauge is the Landau gauge, which can be formulated as solving the following 
conditions in the continuum, 

d^A^ = , (5) 

or as a functional minimisation problem |24j of the Landau gauge fixing functional 

FA[g] =11 A<^ f = - [ d'xTT{{gA,g^f) , (6) 



whose first derivatives with respect to gauge parameters g are shown to give Eq. ([S]) . The Hessian 
matrix Mpp of this functional is the FP operator. 

The functional minimisation approach carries over straightforwardly to the lattice formulation. 
The standard choice of the lattice Landau gauge-fixing functional, in the following called the 
naive or direct lattice Landau gauge functional, to be minimised with respect to the correspond- 
ing gauge transformations gt, is 

Fuig) = E(l - ]^ ReTr5jc/,,^g,+^), (7) 

for SU(A^c) groups. Taking fi{g) = g = for each lattice site i gives the lattice diver- 
gence of the lattice gauge fields and in the naive continuum limit recovers the Landau gauge 
dfiAfj^ — 0. The corresponding Mpp is the Hessian matrix of Fu{g), with respect to the gauge 
transformations. 

Neuberger showed [9l [10] that when all the stationary points, not only minima, of the naive 
functional are taken into account, Zgf turns out to be zero and the expectation value of a 
gauge-fixed variable then is 0/0. Schaden [25] interpreted the problem in terms of Morse theory 
and showed that Zqf calculates the Euler character x oi the group manifold G at each site of 
the lattice. That is, for a lattice with N lattice sites, 

Zgf = Esign(det Mppig)) = {x{G)f, (8) 



where the sum runs over ah the Gribov copies. The group manifold G for compact U(l) is a 
circle S'^ and for SU(A^c) it is 5'^ x 5^ x • • • x 5'^^="^. This result is derived from the Poincare- 
Hopf theorem which asserts that x(M[) is equal to the sum of signs of Hessian determinants 
at all critical points of a non-degenerate height function which is a compact, differentiable and 
orientable function from the manifold M to R. In the case of lattice Landau gauge fixing we can 
immediately identify the gauge-fixing functional ([7]) as a height function, Gribov copies as the 
critical points and Mpp as the corresponding Hessian matrix. 

For compact U(l), G = S^ . As Eq. ([7]) can be viewed as a height function from 5^ x 5^ • • • x S"^ 
to R, and since x{S^) = Oj "^^ have Zqf = 0. As we will explicitly see in section |3l antiperiodic 
boundary conditions in this case fix the global gauge freedom, so the corresponding manifold is 
{S^)^ . For periodic boundary conditions the corresponding manifold is {S^)^~^ . Thus, for any 
boundary conditions, we have Zqp ~ 0. In fact, for any SU(A^c), and hence for the standard 
model group, the group manifolds are odd-dimensional spheres for which x is zero. 
Following this interpretation, for SU(2) gauge theory, Schaden proposed to construct a BRST 
formulation only for the coset space SU(2)/U(1) for which x 7^ (the residual compact U(l) 
symmetry was left unaddressed) . This procedure can be generalised to fix the gauge of an SU( A'^c) 
lattice gauge theory to the maximal Abelian subgroup (11(1))^'="^, since x(SU(A''c)/U(l)) y^ as 
well. This indicates that the Neuberger 0/0 problem for an SU(A^c) lattice gauge theory actually 
lies in (U(l))^'=~^, and hence can be evaded if the problem for compact U(l) is evaded. This is 
why we concentrate our study on the compact C/(l) case from now on. 

We should emphasise here that there are no Gribov copies in continuum Quantum Electrody- 
namics (QED) which is a U(l) gauge theory. More specifically, on the lattice, the compactness 
of the gauge group introduces Gribov copies for the compact U(l) theory. Thus, Gribov copies 
on the lattice are purely lattice artefacts and so is the Neuberger 0/0 problem. As argued above, 
we are interested in studying Gribov copies of compact U(l) theory because the Neuberger 0/0 
problem apparently lies in the residual compact U(l) subgroup. Having said that, we should also 
note that compact QED on the lattice may serve as a prototype for theories based on compact 
groups such as SU(A''c), and hence studying compact QED on the lattice is important in its own 
right. Compact QED in four dimensions exhibits two phases: a Coulomb phase, with a massless 
photon, and a confined phase, which, although unphysical, shares many qualitative features with 
QCD. 

To avoid the Neuberger 0/0 problem, one may modify the gauge fixing condition in some way. 
However, any such modified lattice gauge-fixing should satisfy the following conditions: 

1. The corresponding n\U] should be orbit-independent. If this is the case, then n\U] cancels 
out in Eq. ^ , otherwise it must be computed for each configuration, which is usually not 
feasible as it involves finding all solutions to fi{U^) = 0, for i running over all the lattice 
sites. 

2. It should be possible to efficiently implement the corresponding gauge- fixing numerically. 

3. The additional gauge-fixing terms should not destroy the theory, e.g., the gauge-fixed 
action should be renormalisable. 

Below we discuss a few alternative gauge fixing procedures: 



Minimal lattice Landau gauge: Here, instead of taking all the stationary points of the gauge- 
fixing functional, one only considers the space of minima, called the first Gribov region. Since 
Mpp is positive definite for the minima, by definition, there is no cancellation among the signs of 
determinants of Mpp (i.e., n[U] is just the total number of local and global minima), and hence 
no Neuberger 0/0 is present there. This is also advantageous numerically, since finding minima is 
easier than finding general stationary points. It has also been shown [26| that this restriction can 
be written in terms of a renormalisable action with auxiliary fields (for a review of this approach, 
see [27j). However, a crucial point is whether the corresponding n[U] is orbit-independent. If it is 
orbit-dependent then each orbit comes with a different n[U] and the functional integral becomes 
a function of n[U] and cumbersome to deal with. In the one-dimensional case for compact U(l), 
it has already been shown that the corresponding n[U] is orbit-dependent [551 [55]. In the present 
paper, one of our goals is to verify this in the two-dimensional case. 

Absolute lattice Landau gauge: In this gauge, one further restricts the gauge-fields to 
the space of global minima, called the fundamental modular region (FMR) . Here, in addition to 
evading the Neuberger 0/0 problem by avoiding the cancellation among the determinants of Mpp 
(i.e., n[U] is nothing but the number of global minima), we also expect the corresponding n[U] 
to be orbit-independent, and equal to 1 in the general case. It is also anticipated that the set of 
configurations with degenerate global minima is a set of measure zero which forms the boundary 
of the FMR. In other words, there are no Gribov copies inside the FMR [SUl 131] • In the one- 
dimensional compact U(l) case, this was also verified to be true [28l[29]. In the current paper, we 
want to study this issue in the two-dimensional case. Having said that, we emphasise that finding 
the global minimum of such functions (which corresponds to spin glass model Hamiltonians) is 
known to be a very difficult task and in most cases it is an NP hard problem. Thus, in realistic 
cases, we can not expect to find the global minimum using conventional numerical minimisation 
methods, and the best one may do is to generate a number of minima and choose the 'best 
minimum' among these as an approximation to the global minimum (see eg. [191 132)). In 
addition to this, the absolute lattice Landau gauge can not be stated in terms of algebraic 
conditions, making it difficult to impose using the standard FP procedure. In this paper, we 
provide evidence that n[U] for this gauge is indeed orbit-independent. 

Stereographic lattice Landau gauge: A modification of the group manifold of compact 
U(l), i.e., a circle 5*^, via stereographic projection at each lattice site was proposed and studied 
in Refs. [T31[TSJ[23- For such a stereographically projected manifold the corresponding x is non- 
zero and thus the Neuberger 0/0 problem is completely avoided. Applying the same technique to 
the maximal Abelian subgroup (U(l))^'=~^, the generalisation to SU(A''c) lattice gauge theories 
is possible when the odd-dimensional spheres 5'^'^+'^, fc = 1, . . . ,iVc — 1, of its parameter space 
are stereographically projected to the real projective space RP{2k). It was also shown using 
topological arguments that the number of Gribov copies is exponentially suppressed for the 
stereographic compared to the naive gauge for compact U(l), and that the corresponding n[U] 
is orbit-independent, at least for the one-dimensional lattice case. It can be shown that the 
corresponding FP operator in this case is generically positive (semi-) definite (i.e., all stationary 
points are minima in this case, and there are no saddle points nor maxima) and hence n[U] 
is nothing but the total number of local and global minima. The stereographic lattice Landau 
gauge is a promising alternative from the lattice BRST symmetry point of view since it fulfils 
all the above mentioned practical requirements, except that the orbit-independence is yet to be 
confirmed for lattices in more than one dimension. 



Interestingly, in lattice formulations of supersymmetric Yang-Mills theories, non-compact pa- 
rameterisations (similar to the stereographic projection) of the gauge fields are used |33j . The 
non-compact parametrisation, unlike the compact (group based) parametrisation, apparently 
evades the sign problem in the lattice versions of these supersymmetric theories [341 135] . Re- 
cently, a deeper and direct connection between the sign problem in lattice supersymmetric 
theories and the Neuberger 0/0 problem has been established [35]: essentially, the complete 
action of A/" = 2 supersymmetric Yang- Mills theories in two dimensions can be shown to be a 
gauge-fixing action, analogous to the Faddeev-Popov procedure, with the symmetry now be- 
ing a topological gauge symmetry. Thus, the corresponding partition function is nothing but 
Zgf for the corresponding theory and hence the Neuberger 0/0 problem follows for the com- 
pact parametrisation. For the non-compact parametrisation, the Neuberger 0/0 problem can be 
avoided due to topological arguments. 

Other approaches and efforts: An alternative to fixing the gauge completely would be to 
average over Gribov copies, as proposed in Refs. [371 [SHI 131], with a weight proportional to 
the exp{— PcpFlg]) , where F[g] is the gauge fixing functional of ([5]) or ([7]). This reformulation 
evades the Neuberger problem. In the original formulation, the average was over the whole 
gauge orbit, and the (absolute) Landau gauge would only be reproduced in the limit /3 — >■ oo. 
It has recently been shown AO that if the average is taken only over configurations satisfying 
the Landau gauge condition, the resulting action is renormalisable and possesses a BRST-like 
symmetry. This class of gauges has been studied numerically [41j : however, it is expensive as 
it amounts to performing a second Monte Carlo integral over Gribov copies on top of the usual 
lattice Monte Carlo sampling. 

An orthogonal approach was recently proposed [42], in which the Gribov copies are treated as 
a residual, nonperturbative gauge degree of freedom which (at least in the continuum, infinite 
volume limit) can be fixed by imposing additional conditions on the correlation functions. For 
example, it is suggested that one may define a family of gauges, so-called Landau- -B gauges, by 
defining a target value B for the renormalised ghost propagator at a specific momentum and 
choosing the Gribov copy that gives a result closest to this value. A particular choice is the 
max-i? gauge which selects the Gribov copy with the largest value for B. Whether or to what 
extent it is possible to impose such a condition is at present not clear, and little or nothing is as 
yet known about the relation between the B parameter and the Landau gauge fixing functional. 
Another way was put forward in Refs. [i31 H^ H5] using the ghost/anti-ghost symmetric Curci- 
Ferrari gauges . There the argument used was that the Neuberger 0/0 problem could be extended 
to include such non-linear gauges with their extended double-BRST symmetry despite their 
quartic ghost self-interactions, which allow the introduction of a mass term for ghosts. Such a 
Curci- Ferrari mass would break the nilpotency of the BRST/anti-BRST charges which is known 
to result in a loss of unitarity; however, this mass also serves to regulate the Neuberger zeroes 
in a lattice formulation and expectation values of observables can then be meaningfully defined 
in the limit m — ?• via rHospital's rule. 

On the continuum side, an explicit counting of Gribov copies was carried out in J20U21] . There, 
the SU(2) case was considered, but restricted to static spherically symmetric configurations only 
(in which case Landau and Coulomb gauge conditions are identical). 

This is not intended to be an exhaustive review of previous attempts to address this subject. 
For a pedagogical and thorough review of previous work on this topic, we refer to the recent 
review 1461. 



3 Lattice Landau gauge for compact U(l) 

Since the group manifold of compact U(l) is topofogically a circle S^, we can write the link 
variables and gauge transformations in terms of angles 4'i^^,0i G (— tt, tt] mod 27r, as Ui^^ = 
e"^*'^and gi = e*^% respectively. Thus, the naive gauge fixing functional Eq. ([T]) is reduced to 

^<^(^) = ^ E (1 - CO^('^^M + e^+f. - 0^)) ^ :^ E(l ~ COS0,?^), (9) 

where we have defined (f>^ :— 4>i,^ + di+fi^0i. Note that we have also introduced a normalisation 
factor l/V, where V is the lattice volume. The stereographic gauge fixing functional becomes 

2 
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F;{e)^-^Y.^n{cos{<j>tj2)). (10) 



A given random set of (j^i^fj, is called a random orbit or a hot configuration. The special case when 
all (/)i,^ are zero is called the trivial orbit, or cold configuration. The gauge-fixing conditions are, 
respectively, 

d 

M(^) = -J] (sin <^ -sin C/,, J -0, (11) 

d 
ftiO) - - E (tan«^/2) - tan(CA,p/2)) - (12) 

for all lattice sites i. The FP operator for the two gauges are 



iMFp)i,j = E ( " CO^(t>l^'>^+^^,J + (COS-^-,^ + C0S<^f_^_^)J,j - COS( 



A 



h-i,,^o,^fiJ) ' (13) 

(M'^ph. = E ( - '^^" 4^^^+A,. + (sec^ ^ + sec^ ^)5,, - sec^ %^5._,,) . (14) 
The boundary conditions are given by 

Oi+Nf, = {~l)'"0„ (t>i+Nil.t, ^ (-1)V-*,M> (15) 

where TV is the total number of lattice sites in the /i-direction. We have fc = for periodic 
boundary conditions (PBC) and A: = 1 for anti-periodic boundary conditions (APBC). With 
PBC there is a global degree of freedom leading to a one-parameter family of solutions, as all 
the equations are unchanged under Oi ^- Oi+ -d, Vi where -d is an arbitrary constant angle. We 
remove this degree of freedom by fixing one of the variables to be zero. In practice, we set 

d{N.,....N) = 0. 

From now on we concentrate only on these gauge-fixing functionals and let the {0i,^} take 
random values independent of the action. This corresponds to the strong coupling limit (3 — 0. 
Note that this is sufficient to answer the question of whether or not n[U] is orbit dependent, as 
every gauge orbit has a non- vanishing weight for any finite /3. 

It is worth noting that the functional ([9]) is identical to the hamiltonian of the random phase 
XY model, and that the functional for the trivial orbit is identical to the hamiltonian of the 
classical XY model in statistical physics [HI [21] . 



3.1 Analytical results in 1 dimension 

Here we list the available analytical results in one dimensiono 

1. For the naive functional with APBC I2S1II1], the minima for any orbit are 0^ = or tt for 
all i — I, . . . , N, the number of minima for any orbit is 2, and the number of stationary 
points is 2^. 

2. For the naive functional with PBC [551 [55] the stationary points for any orbit and any odd 
number of sites iV, are given by 

0,^ = (-!)«' </)?^+q,7r mod27r, q, e {0,1} , (16) 

with (j) := j^ X]i=i 4'i ■ In fact, since we are interested in solutions of (/)f only modulo 
multiples of 2tt, it is sufficient to consider 
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The number of stationary points is X]7=o 1-^ ~ 2j|( 7 j, i.e., it increases exponentially 
with N . The minima occur when qi = 0, for alH = 1, . . . , A^ — 1, and cos 4)% > 0. Thus, the 
number of minima is bounded by N . The number of Gribov copies for the compact C/(l) 
theory in one dimension thus increases exponentially, but the number of Gribov copies in 
the first Gribov region increases only linearly. In [5S], it was also shown that if = + ^)7r 
with j G N, the FP operator is singular, i.e., for these orbits Gribov horizons do exist. F 
evaluated at these stationary points is 

N 

fe=i 

The minima of this function are when both conditions (1) q^ =0 for all fc = 1, . . . , A^ — 1, 
and (2) cos(/)^ > with i/)^ given as Eqs. ([T6|) and ([T7|. are satisfied. Thus, the minima 
are given by 

It follows from (1201) that function values at the minima accumulate near zero as iV — ^ c». 
The proof of this is given in Appendix [C| 

3. We show in Appendix [K\ that for the stereographic gauge with APBC, the FP operator is 
generically positive definite and hence the function F^iO) has only minima. In fact, there 
is only one minimum for any orbit in this case, (j)f = mod 2tt, i = 1, . . . , N. 



^We note that some analytical results for both U(l) and SU(2) in 1 dimension were already found in I47| . 



4. For the stereographic gauge with PBC, again, the FP operator is generically positive 
definite as shown in Appendix \K\ There are N minima for any orbit in this case |28[|48] . 

.a -7 ^TT 



-r mod 2^, r = 0,...,^-l, ^ - ^ I] "^^ ' (^1) 

4=1 

i.e., the number of minima is orbit- independent and increases hnearly with N . 



3.2 Some analytical results in 2 dimensions 

1. In any dimension, for the trivial orbit, there is only one global minimum with PBC, and 
two global minima with APBC, for both the naive and stereographic gauge. 

Proof: Consider © defined on a lattice of size N"^. Let the components of 9 be given 
by {^(ii,...,ij)} where each 9(^i^i^-^ £ (~^i ""li ii, ■ ■ ■ jid runs from 1, . . . , iV, and d is the 
dimension of the lattice. It is clear that the minimum value the naive Landau gauge-fixing 
functional, given in Eq.([9]), can attain is 0. Since it is a no n- negative function, when F = 0, 
the only possible configuration is when each term in F is zero, i.e., 



7(ii,...,jrf)+/i - f(ii 



= V^l,...,^dG{l,...,7V},/ie{l,...d}, (22) 



which is the global minima. The stereographic Landau gauge-fixing functional is also a 
non-negative function and hence all the above arguments apply. Then any global minimum 
has to satisfy (|22|) and so with 9 defined as above we must have 0(ii,...,i^)+/i = ^(ii,...,^^) 
mod 27r. Therefore, all 6'(ij^ i^) can be parameterised by only one of them, say, ^(i....^i) = 
^(ii,...,id)+A' fo^ ^11 ii,-.-,id and /i. 

Now, for the APBC case, since 9(^ii]y^i-j — —9(^ii-j but also from the above arguments 
^(i,...,i.Af+i) = ^(i....,i)niod27r. So we have 29(^ii) — mod 27r. Thus either 6'(i i) = 
or ^(i,...,i) = TT. So all global minima of © with APBC are lattices where all of the 
^-variables are zero or all 6'-variables are tt. With PBC, we must fix one element on the 
lattice and require it to be zero. But as all elements are equal, this means that all the 
elements must be zero and thus the only global minimum of ([9]) with PBC is a lattice with 
all 0- variables set to zero. 

2. For a generic orbit for compact U(l), the number of stationary points for the naive func- 
tional is > 2^ for APBC and > 2^^^ for PBC. For the stereographic gauge, this lower 
bound is 1. For generic SU(A^c) with PBC, the number of stationary points is lower- 
bounded by 2(^<^~^)(^~^) (See Appendix IB] for a proof). 

3. If is a minimum of the naive or stereographic lattice Landau gauge with APBC, then 
one can construct another minimum by adding n to every 9i. 

This is easy to see in one dimension: let 6*^ -^ 9i + n, for i — 1, . . . , TV — 1, then (pf ^> (pf. For 
i = N, this transformation leaves 0^ = (J)n — 9i—9n —^ 4>n ~ 9i — 9n — 2tt = 0^ mod 2tt. 
Due to the periodicity of the trigonometric functions there is no effect of this transformation 
in either the gauge-fixing equations themselves or in the hessian. It is straightforward to 
extend the same argument to any dimension. 



10 



4. For the trivial orbit in d dimensions, for a symmetric N'^ lattice, ii 9 = {^^1,...,^^} is a 
minimum, then 9 = iG9 is also a minimum, where G is any lattice symmetry operation, 
including lattice translations, rotations and reflections (or axis permutations). This follows 
straightforwardly from the symmetries of the defining equations. Specifically, for example 
in two dimensions, ii 9 — {9ij} is a minimum of the naive or stereographic gauge fixing 
functional, then 9'^ = {9ji} is also a minimum, as the following argument shows. 

Let 6* be a solution of Eq.®, i.e., it satisfies the following equations for all i,j ~ I, . . . ,N, 



OF 









ly J- 




= sm{9i,j - 


-0,- 


-1,.) - 


- sin(6'i+i J - 


'^i,i) 




d9i 


J 












+ sm{9ij 


-0, 


..-l) 


- sin(6'j_j+i 


-e^,) 


relabel i 


-^ 


JJ 


-^ 





— am{9j^i - 


^e,- 


-m)- 


SlTL(9.j + l,i - 


- e,,) 












+ sm{9j^. 


-e, 


.-l) - 


- sin(6'j,j+i 


- 0,,) 












OF 
09,., 


= 0. 









(23) 

(24) 



Thus, 0^ is a stationary point of ©■ To prove that O'^ is a minimum, we note that the 
Hessian (TTS)) for the naive and similarly for the stereographic gauge in the case of the 
trivial orbit are both symmetric under the interchange of coordinates, and therefore the 
hessian evaluated at 9 has the same eigenvalues as the hessian evaluated at 9"^ . Thus, if 9 
is a minimum, then 9^^ is also a minimum. 

It is easy to see that this permutation symmetry holds for any dimension, and works for 
both naive and stereographic gauge fixing functionals and for stationary points in general, 
for the trivial orbit. Since there are AN'^'d). elements of the symmetry group for a N'^ lattice 
[N'^ lattice translations, d\ permutations, plus reflection and sign change), this means that 
if is a stationary point of a d-dimensional functional then the function value at 9 will be 
47V'^d!-fold degenerate unless 9 maps onto itself under a subset of the symmetry operations. 
For example, in d = 2 the only permutation of spatial coordinates is a transposition, giving 
a degeneracy of 8A^^, but is reduced by a factor two for any skew-symmetric stationary 
point, 6'"^ = —9. The function value of zero is non-degenerate, ss. 9 — maps onto itself 
under all symmetry operations. 

If we choose to impose periodic boundary conditions, then we essentially choose the fixed 
site to have lattice coordinate in our case {N, . . . , iV), which is fixed under any permutation 
and sign flip, but not under lattice rotations or reflections. This means that the generic 
degeneracy is reduced to 2d\. 

4 Numerical polynomial homotopy continuation method 

In general, systems of non-linear equations are extremely difficult to solve. However, if the 
non-linearity in the system is polynomial-like, then the situation is enhanced due to the recently 
developed algebraic geometry methods. In particular, we will use the so-called numerical poly- 
nomial homotopy continuation (NPHC) method [49] to find all the solutions of the gauge-fixing 
equations. This method was introduced in particle physics and statistical mechanics first in j28j 
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and applied in[Sni[iIl[il[Sl[Sl[i51[SSl[571[SS]. Below, first we show that the problem of solving 
the extremising equations in terms of the 0- variables can be transformed into that of solving 
a system of multivariate polynomial equations. Then we describe the numerical polynomial 
homotopy continuation method, which can be used to find all the solutions of a given system of 
polynomial equations numerically. Finally, we give our results for the problem at hand. 
To convert the naive gauge-fixing equations, we can first use trigonometric identities to rewrite 
(|TT|) as 

fi{c,s) = ^ {c^{ci+f,sm(f)i^^-c^^f,sin(f)i^p,^fJ_ + Si^f,cos(|)i^f,^fJ_ + Sl+f,cos(|)i,fJ.) 

M (25) 

+ s,{si+f, sin (j)i^f^ - Si- f, sin (j),^f,^^ - Ci+/iC0S(/)j,p - Ci_/i cos^i-/^,^)), 

where we have written Si :— sin 9i and Ci :— cos 9i. This is merely a change of notation. However, 
we can now add additional equations to the system for each site i, namely, 

g,{c,.s)^sf + cf-1^0, (26) 

The combined system of all /i(c, s) and gi{c, s) is not just a change of notation: all the Ci and 
Si are now algebraic variables and the equations are multivariate polynomial equations, i.e., the 
fact that Ci and Si are originally sin 9i and cos 6i is taken care of by the constraint equations 
(PS]) . In general, for a lattice with N lattice sites, we have in total 2 A'' polynomial equations and 
2N variables. 

To convert the gauge-fixing equations arising for the stereographic gauge, we first simply expand 
(jl2[) using the trigonometric identity 

X + y + z sin x + cos z sin y + cos y sin z 

tan ^^ ; ; — , 

2 cos X + cos y cos z — sin y sin z 

to obtain 

rs( ] — \^ ( ^^'^ '^'".M^' ~ "^os (j)i,fj,Si + Si+f, sin (j)t_f^^f,Ci-fi - cos cp^-fj^f^Si-f, + Si . 

' ' -^ ^sin0i,^Si -|-cos(/)i,^Ci -t-Ci+/i sin0i_/i,^Si_/i -I- cos0i_/i^pCi_/i -I- Ci ' ^ -' 

which can again be translated into the same polynomial form as above. Here, the difference is 
that the above equations are not in the 'polynomial form' due to the denominator. We can clear 
the denominators out by multiplying them with the numerators appropriately and assuming 
that none of the denominators are zero (such solutions can be sorted and thrown out once all 
the solutions are obtained). 

For example, for the trivial orbit on a one-dimensional lattice with A^ = 3 and APBC, Eq. (P7)) 
simplifies to 

S2 - Sl _ Si + S3 _ S3 - S2 _ S2 ~ Si _ -Sj - S3 _ S3 - S2 _ , , 

Cl + C2 Ci + C3 C2 + C3 Ci + C2 Ci + C3 C2+ C3 

After clearing out the denominators, the equations become 

-2cisi - C2S1 - C3S1 + C1S2 + C3S2 - C1S3 - C2S3 = , 

C2S1 + C3S1 - C1S2 - 2C2S2 - C3S2 + C1S3 + C2S3 = , 

-C2S1 - C3S1 -I- C1S2 -I- C3S2 - C1S3 - C2S3 - 2C3S3 = , 



1 - y(ci -I- C2)(C1 + C3)(C2 -t- C3) = , 
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where the last equation is added to ensure that the denominators are never zero, and y is an 
additional variable. Thus, the final set of the equations is in polynomial form. One can then 
solve this system using the NPHC method. 

4.1 The method 

Let us consider a system of multivariate polynomial equations, say P{x) = 0, where P{x) = 
{pi{x), . . . jPrnix))"^ and x = (xi, . . . ,Xm)'^, which is known to have isolated solutions, e.g., the 
above mentioned gauge fixing equations after eliminating the global gauge freedom. Now, the 
Classical Bezout Theorem asserts that for a system of m polynomial equations in m variables, 
for generic values of coefficients, the maximum number of solutions in C™ is Y\T=i '^«' "where di 
is the degree of the ith polynomial. This bound, the classical Bezout hound (CBB), is exact for 
generic values (i.e., roughly speaking, non-zero random values) of coefficients, e.g., for the one- 
dimensional naive (or minimal) gauge fixing equations with N of lattice sites and with APBC, 
this number is 2^^ (because there are 2N polynomials each of which is a degree 2 polynomial). 
The genericity is well-defined and the interested reader is referred to Ref. [49j [59] for details. 

Based on the CBB, a homotopy can be constructed as 

H{x,t)=-i{l-t)Q{x)+tP{x), (30) 

where 7 is a random complex number. Q{x) ~ {qi{x), . . . , qm{x))'^ is a system of polynomial 
equations with the following properties: 

1. the solutions of Q{x) — H{x,0) = are known or can be easily obtained. Q{x) is called 
the start system and the solutions are called the start solutions, 

2. the number of solutions of Q{x) = H{x, 0) = is equal to the CBB for P{x) ~ 0, 

3. the solution set of H{x, i) = for < i < 1 consists of a finite number of smooth paths, 
called homotopy paths, each parameterised by t G [0, 1), and 

4. every isolated solution of H{x, 1) = P{x) = can be reached by some path originating at 
a solution of H{x, 0) — Q{x) = 0. 

The start system Q{x) ~ can for example be taken to be 

Q{x)^\ : =0, (31) 






where di is the degree of the i*'' polynomial of the original system P{x) = 0. Eq. (PT|) is easy 
to solve and guarantees that the total number of start solutions is rii=i'^j' ^^^ ^^ which are 
non-singular. 

One can then track all the paths corresponding to each solution of Q{x) = from i = to 
t = 1 and reach P{x) = = P[{x, 1). By implementing an efficient path tracker algorithm, e.g., 
Euler predictor and Newton corrector methods, all isolated solutions of a system of multivariate 
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polynomials system can be obtained. The complex random number 7 is crucial here: it has 
been shown [49] that for a randomly chosen 7, there are no singularities (i.e., paths do not cross 
each other) for t € [0, 1). This ensures that in the end we get all the solutions. In this respect, 
the NPHC method has a great advantage over all other known methods for finding stationary 
points or minima. 

There are several sophisticated numerical packages well-equipped with path trackers such as 
Bertini[60], PHCpack [61], PHoM [62] and HOM4PS2 [63l[59]. They ah are available as freewares 
from the respective research groups. We mainly use Bertini, HOM4PS2 and PHCpack to get the 
results in this paper: for each of the systems, we use at least two of the packages to cross-check 
the results. 

4.2 Results 

Before proceeding to the results, it should be noted that a solution here means a set of values 
of Si and Ci (for the naive or minimal gauge) or t^'s (for the stereographic gauge) satisfies all of 
the equations with tolerances lO^^'^. All the solutions come with real and imaginary parts. A 
solution is a real solution if the imaginary part of each of the variables is less than or equal to the 
tolerance 10"^ (below which the number of real solutions does not change, i.e., it is robust for 
all the cases we consider in this discussion). The original trigonometric equations are satisfied 
with tolerance 10^^° after the si and q (or ti) are transformed back to 9i. All these solutions 
can be further refined to an arbitrary precision. 

We present the results for the two-dimensional naive functional by classifying the obtained 
solutions in terms of the number of positive and negative eigenvalues of the corresponding 
Hessian matrix, or the FP operator, because then we can use the Neuberger zero as a necessary 
condition for having all the solutions. Before proceeding to the two-dimensional case, we note 
that we have reproduced the known analytical results in the one-dimensional case for both APBC 
and PBC naive and stereographic gauge cases using the NPHC method, up to TV = 25 for various 
random orbits. 

We now explore the simplest non-trivial case in higher dimensional lattices which is the naive 
functional on a 3 x 3 lattice with the trivial orbit (TO) and 10 random orbits (ROs) (i.e., 
randomly chosen (j>i^^ G (— 7r,7r]) with PBC and APBC. For these cases, the CBB values are 
65536 and 262144, respectively. 

The results of our runs are summarised in Tables [1] and [2] for APBC and PBC, respectively. We 
find that: 

1. For the trivial orbit with APBC, there are a total of 3768 real solutions, i.e., Gribov copies. 
Out of these, there are 1952 solutions which have zero FP determinant with tolerance 10~®. 
These singular solutions lie exactly on the Gribov horizon. They can be further classified 
in terms of the number of zero eigenvalues of the FP operator at each of the solutions. 
This amounts to classifying the singular solutions of the polynomial system in terms of 
their multiplicities using the so-called deflation singularities technique [63] . The remaining 
1816 real solutions have non-singular Mpp- Similarly, the trivial orbit with PBC has 30 
singular (Gribov horizon) solutions out of a total of 1112 Gribov copies. 

2. There are no singular solutions for any of the random orbits for either PBC or APBC. 
The total number of Gribov copies for each random orbit fluctuates around ^ 2500 for the 
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54 


10 
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44 
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44 
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2672 


10 


60 


208 


460 


598 


598 


460 


208 


60 


10 


R07 


2504 


2504 


8 


46 


190 


426 


582 


582 


426 


190 


46 


8 


R08 


2304 


2304 


6 


48 


152 


362 


584 


584 


362 


152 


48 


6 


R09 


2352 


2352 


6 


50 


182 


384 


554 


554 


384 


182 


50 


6 


ROlO 


2382 


2382 


6 


44 


180 


402 


558 


558 


402 


180 


44 


6 



Table 1: Summary of the Gribov copies of the naive lattice Landau gauge for 3 x 3 lattice, 
for different orbits with APBC. A^tot is the total number of real solutions and Nns the number 
of non-singular solutions. RO denotes random orbit while TO denotes the trivial orbit. The 
number of global minima is 2 for all orbits. 
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Table 2: As Table [l] but for PBC. TO mag field denotes the trivial orbit with the modified 
functional ([32]) with h = 0.01 (see main text for explanation). The number of global minima is 
1 for all orbits. 
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APBC case and around ~ 500 for the PBC case. We see that the total number of Gribov 
copies is an orbit-dependent quantity in both cases. 

3. All the Gribov copies (excluding the Gribov horizons) can be classified by the number of 
negative eigenvalues, as shown in Tables [1] and [H In the APBC case, one can see a perfect 
symmetry among the Gribov copies in the number of solutions classified according to the 
negative eigenvalues of the FP operators evaluated at these solutions. This symmetry yields 
a perfect cancellation of the signs of the FP determinants, giving rise to the Neuberger 
zero. For the PBC case, there is no such manifest symmetry. However, we can easily check 
that the sum of signs of the FP determinants is still zero for all orbits, even though the 
total number of Gribov copies is an orbit-dependent quantity. 

4. The solutions with no negative eigenvalues are the minima. We see that for both types 
of boundary conditions, the total number of (local and global) minima is orbit-dependent. 
The global minima can be determined by computing the functional F with appropriate 
boundary conditions for each of these minima and identifying those giving the lowest 
value for F. We find that the number of global minima is orbit-independent. It should 
be emphasised that though there are two global minima for the APBC case, both of them 
are trivial copies of each other. 

5. Strictly speaking, because there are singular solutions for the trivial orbit with both APBC 
and PBC, Morse theory does not directly apply hereO Hence one should not expect that 
the sum of the signs of the FP determinants would give the Euler characteristic in such a 
case, and this can indeed be seen for the trivial orbit with PBC in Table[2] Instead, we have 
to work with Morse-Bott theory, a discussion of which is beyond the scope of this article. 
Note that Schaden's equivariant construction takes this fact into account. However, a 
conventional way to remove such singular solutions is to add an external magnetic field 
term in the original height function, which in our case amounts to adding an external 
magnetic field term to Eq. ([9]), i.e., 

N 

F^{0) ^ F^^e) = F^{e) + hY^ cose,. (32) 

The new system of equations can again be translated to polynomial form and solved using 
the NPHC method. The resuhs for the trivial orbit with PBC (taking h = 0.01) are 
also reported in Table [U where one can again see that the sum of the signs of the FP 
determinants is zero. 

From our results we can also find the distribution of functional values F, broken down by the 
number of negative eigenvalues (ie, for each Gribov region). The results are shown in Fig. [T] 
for antiperiodic boundary conditions; the results for periodic boundary conditions are similar. 
We see that the functional values cluster around different values for different n, with a peak of 



^It is important to mention that we have observed using the numerical algebraic geometry methods that there 
is a positive dimensional real component of dimension at least one in the solution space of 3 X 3 lattices. This is 
surprising because after eliminating the global gauge degree of freedom, there should only be isolated solutions 
left. Further studying this phenomenon and whether this is of any physical significance is very interesting but 
beyond the scope of the present article. 
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AF = F 



Figure 1: Distribution of the values of the gauge fixing functional relative to the absolute 
minimum for each orbit, AF — F ~ -Fmin, for 10 random orbits in the naive gauge on the 
3x3 lattice with antiperiodic boundary conditions, and for different number n of negative 
eigenvalues. Top: Histogram of the absolute number iVcopy of copies found in each bin; Bottom: 
The distribution separately normalised for each value of n. 

the distribution which is roughly linear in n. However, there is a substantial overlap between 
the distributions for different n, and in particular (because of the much larger number of saddle 
points than of minima) there are many saddle points with the same functional value as that of 
typical minima. 

It would clearly be interesting to find out how this pattern changes for larger lattice sizes, but 
this is beyond the scope of this study. 

The equations for the stereographic gauge are quite difficult to solve: the polynomial form of 
the equations is both dense and of high degree. We have only been able to find all the solutions 
for the trivial orbit on the 3x3 lattice with periodic boundary conditions, where the number of 
Gribov copies was found to be 5256. This is much larger than the 112 copies found in the naive 
gauge, showing that the exponential suppression in 1 dimension does not carry over to higher 
dimensions. 
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5 Numerical minimisation 

So far we have collected or worked out some analytical results for solutions in the one dimensional 
lattices and have used the NPHC method to find all the solutions for the 3x3 lattices. For 
lattices larger than 3 x 3, we have not yet been able to obtain results using the NPHC method. 
To get results for the bigger lattices, we have to rely on the traditional 'brute force' methods. 
Note that for these traditional methods, there is no guarantee of finding all solutions nor all 
minima of a given function. But our strategy is to first study the cases for which we already 
know all the solutions: we reproduce the results for the one dimensional lattices and for the 3x3 
lattice using the traditional method. This gives us a hint on how much computational effort we 
will require for larger lattices. 

In particular, we have used the conjugate gradient algorithm [65) to numerically obtain minima 
and stationary points for both the naive and stereographic gauge. Once the conjugate gradient 
algorithm has converged we check that || VgF |p< ei where ei is sufficiently small to ensure 
that we have actually obtained a minimum. We reduce 9 mod 2tt to make sure the values are 
in the interval (— tt, tt]. We have to set a tolerance, 62, on how close 6i has to be to tt or 
to be considered either, respectively, and apply these changes. Errors occurred when trying to 
determine the total number of minima if this flag was not in place. Also, if \9i + 7r| < £2 then 
di —7- TT, since 9i e (— tt, vr], not [— tt, tt]. 

To obtain "all" minima, we generate a series of random initial guesses for the variables Oi, and 
successively minimise the given functional from each initial guess. We can then count the number 
of unique minima. After increasing the number of samples sufficiently, it is noted that in most 
cases the total number of unique minima that the function converges to stops increasing and 
thus we can say to have completely sampled the solution space of the function. Examples of this 
are shown in Fig. [2l We note that the number of unique minima found is usually far smaller 
than the number of samples used. We call this method the Monte Carlo Conjugate Gradient 
Method (MCCGM). 

5.1 1 dimension 

Since our method does not guarantee that we find all the minima, we first show that we reproduce 
all the minima in the one-dimensional case where rigorous analytical results are available. Indeed, 
we have reproduced all the minima for lattices up to iV = 1000, for both the naive/minimal 
and stereographic gauge with APBC. We also reproduced all the minima with PEC for up to 
A^ = 16 for the naive functional and up to A^ = 7 for stereographic gauge. In Figure [31 we 
plot the number of minima as a function of the number of lattice sites for both gauges with 
both types of boundary conditions With APBC the number of minima is always 2 for both the 
stereographic and naive gauge functional. With PBC, the number of minima is orbit-dependent, 
and correspondingly two possible values are shown for each A^. These results are in complete 
agreement with the analytical results found previously. 

To separate out the global minimum, we have then evaluated the corresponding functional value 
at each of the local minima, and identified the minimum (or minima) with the lowest functional 
value. To remove any uncertainty, we have also used Matlab's in-built routines for Simulated 
Annealing and Genetic Algorithm to find the global minimum for the one-dimensional cases, and 
these global minima always match the one obtained from the conjugate gradient minimisation. 
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Figure 2: Number of distinct minima -/Vmin found as a function of the number of random initial 
guesses A^guess (start vectors for the CG minimisation), for 9 different random orbits in the naive 
gauge on the 5x5 lattice with periodic boundary conditions. Note the logarithmic scale on the 
horizontal axis. 

The global minima from all these methods match up to A^ = 35 for the naive APBC case. 
For lattice sizes larger than this the Matlab minimisation routines became inefficient both in 
terms of memory and computation time. The MCCGM was found to be superior to both of 
these methods as it took less time and memory to find all minima compared to the Simulated 
Annealing and Genetic Algorithm routines. We find that for all the random orbits, the number 
of global minima is the same (one for PBC and two for APBG). 

In figure m we show the distribution of function values at each minimum for 100 different orbits, 
for the one-dimensional naive (minimal) gauge with PBC, for four different values of the lattice 
size A^. We see that for small A^ the function values are almost uniformly distributed between 
and 1, but as A^ grows the function values tend to accumulate near 0. This is in accordance 
with the arguments in Appendix [C] that function values accumulate near zero as A^ — )■ oo. 

For the stereographic gauge, our numerical minimisation approach broke down beyond A'^ = 16 
for PBC, in that we failed to find all the minima even when increasing the number of initial 
guesses to beyond 10^ (for APBC it worked fine up to A^ = 1000). In this case each minimum 
is located in infinitely deep potential wells, and which particular minimum we find is uniquely 
determined by the initial guess. A strategy to select the initial guesses by uniformly sampling 
these potential wells in the PBC case has proved elusive. The same problem affects the 2- 
dimensional stereographic gauge even more strongly, as we shall see in the next section. 
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Figure 3: The number of (local and global) minima A^min vs the number of lattice sites N, for 
naive (left) and stereographic gauge (right) in 1 dimension, with PBC and APBC, for the trivial 
orbit (TO) and random orbits (RO). 



5.2 2 dimensions 

5.2.1 Naive lattice Landau gauge 

Table [3] and Fig. [5] show the number of minima Ny^in for the trivial orbit of the naive functional 
with both boundary conditions. Because of the permutation symmetries noted in Sec. 13.21 each 
minimum will in general be n-fold degenerate, where n = 16N^ for APBC and n — A for 
PBC. For smaller lattice sizes in particular, some of the minima are invariant under a subset 
of these permutations, reducing the level of degeneracy. For periodic boundary conditions, 
we find a higher level of degeneracy than naively expected, suggesting additional symmetries. 
Indeed, studying the individual solutions we find that they consist of permutations of the same 
numbers. We have been able to classify some of these permutations, but have not as yet found 
any systematic pattern that is valid for generic lattice sizes. 

It is clear that the number of distinct minima increases considerably for N > 7, but it is not 
possible to conclude from our results on these small lattices whether it increases polynomially 
or exponentially with N. 

In FigurelHwe plot the function values at each minimum against N for each lattice of size N x N 
up to a lattice of size 12 x 12. This shows that the function values of 2-D minima of the trivial 
orbit appear to decay towards zero, just as in the 1-D case. This is again in accordance with 
the arguments in Appendix O For the lowest non-trivial minimum, which clearly can be seen to 
decrease with N, we see by inspection that this consists of a layering of 1-dimensional minima, 
as described in Appendix \C\ whose function values (at least for the lower-lying ones) can be 
shown to decrease with N. However, we do not have a proof that most minima approach a 
function value of zero. We shall see later that this behaviour does not carry over to random 
orbits in 2 dimensions. 

To study the random orbit case, we have generated iVorb random orbits for each lattice size. For 
each orbit, we have found all the solutions, including the global minimum. The details are given 
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Figure 4: Histogram of function values at each minimum for 100 different orbits, for the one- 
dimensional naive functional with PBC. 

in Table m 

Our first observation concerns the number of solutions (Gribov copies), which, as expected, is 
orbit-dependent and increases with N. Unlike the trivial orbit case, all the minima are distinct 
for the random orbits (or twofold degenerate for antiperiodic boundary conditions). The number 
of minima for the 3x3 lattice agrees with that found using the NPHC method in section|4l Note 
that because of the twofold degeneracy for APBC we have divided the total number of solutions 
by 2 to factor this out. As we can see in Fig. [7] the number of solutions increases exponentially 
with N, with roughly the same rate of increase for periodic and antiperiodic boundary conditions. 
We also found that there is always 1 global minimum for PBC and 2 degenerate global minima 
for APBC. For a more detailed view, figure |8] shows histograms of the number of solutions for 
the four smallest lattice sizes. We do not observe any significant change in the shape of the 
distribution as N increases, at least for these small lattices. 

We now turn to the values of the gauge fixing functional F for the different lattice sizes. Two 
issues are of particular interest: firstly, how does the absolute minimum change as a function 
of N; and secondly, how does the spread of function values change as A^ changes? If the first 
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Figure 5: Plot of the lattice size TV vs number of (global and local) minima of the naive functional 
in 2D, with PBC and APBC, for the trivial orbit. 

Gribov region becomes equivalent to the fundamental modular region in the large volume limit, 
as conjectured, then we might expect the functional values to accumulate closer to the absolute 
minimum for larger volumes. We therefore determine the absolute minimum as well as the 
median value of the functional at the minimum for each random orbit. The results are shown 
in Table m 

Unlike the case for the trivial orbit, we find no evidence that the absolute minimum -Fmin 
decreases towards zero as N is increased. Instead, i^min appears to converge to a value around 
0.4. This is not unsurprising, as we are working with random gauge configurations (corresponding 
to the strong coupling limit), which are not smooth and should therefore not be expected to yield 
smooth gauge fixed configurations characterised by a small value for the gauge fixing functional. 

We also find no clear evidence that the typical (or median) functional value at the minima ap- 
proaches the absolute minimum as TV increases. If anything, the indication is that the opposite 
is the case for periodic boundary conditions, while for antiperiodic boundary conditions the me- 
dian value may be getting closer to the absolute minimum. This is borne out by the distribution 
of functional values relative to the absolute minimum, shown in Fig. |9l For both boundary 
conditions the distribution goes from being fairly flat to more sharply peaked as N increases, 
but for APBC the peak appears to shift towards zero, while for PBC the peak position stays 
constant or increases slightly. 

As an aside, we observe that the lower minima are always the more likely to be found by the 
minimisation algorithm, with the global minimum among the first to be found. If instead of 
counting each Gribov copy equally, we weight them with how frequently they are found by the 
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Lattice Size N 


3 


4 


5 


6 


7 


8 


9 


10 


Nmin (APBC) 


2 


66 


202 


146 


1570 


7170 


24626 




(iV„,i„ - 2)/N^ 





4 


8 


4 


32 


112 


304 




-^distinct 


1 


2 


2 


2 


4 


13 


20 




^min (PBC) 


1 


1 


5 


25 


9 


671 


>4400 


> 12000 


iVdistinct (PBC) 


1 


1 


2 


4 


3 


9 


14 


>39 



Table 3: Details about the total number of minima A^min and distinct (non-degenerate) minima 
A'distinct for the naive functional, trivial orbit, on N x N lattices. 



conjugate gradient search, we find a distribution which is peaked at zero and appears to become 
more strongly peaked as N is increased. From a practical point of view this suggests that 
numerical minimisation methods may have a reasonable chance of approaching the fundamental 
modular region even for larger lattices, without having to explore more than a small fraction of 
the Gribov copies within the first Gribov region. 

5.3 Stereographical lattice Landau gauge 

We have also attempted to use the MCCGM to find the minima of the 2-D stereographic gauge 
fixing functional. In this case, we have been restricted to a 3 x 3 lattice with APBC, and even 
there we have not succeeded in finding all the minima, except for the trivial orbit. However, we 
clearly see that the number of Gribov copies (minima) for the stereographic gauge is larger than 
the total number of Gribov copies (all stationary points) for the naive lattice Landau gauge. 
Moreover, the results we have obtained suggest that n[U] for the stereographic gauge is orbit 
dependent. As an illustration of this. Fig. [10] shows the number of distinct minima found for 
three random orbits, as a function of the number of initial guesses. We find that the number of 
Gribov copies is greater than 4500 for both random orbits considered, and 7466 for the trivial 
orbit. We may compare this with the total number of Gribov copies for the naive gauge shown in 
Table[TJ which ranges between 2300 and 2700 for the 10 orbits considered there. This shows that 
the exponential suppression of the number of Gribov copies observed in 1 dimension does not 
hold in higher dimensions, and the stereographic gauge is thus less advantageous from this point 
of view than the naive gauge. Furthermore, the number of minima found for the two random 
orbits and the trivial do not appear to converge on a single number, although this cannot yet be 
completely ruled out, as the convergence for the random orbits seems only logarithmic in A^gucss- 
Similar results were obtained using the NPHC method for 2 x 2 lattices with periodic boundary 
conditions. 
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Figure 6: Function values at the minima for the 2-dimensional naive functional, trivial orbit, as 
a function of the lattice size N. Note that the global minima at F = are not shown here. 



6 Discussion and conclusions 

In this paper, we have studied three alternative definitions of the lattice Landau gauge: minimal 
lattice Landau gauge, absolute lattice Landau gauge, and the stereographically modified lattice 
Landau gauge. We have focused on the compact U(l) case only for this paper, with (3 — 0. 

Our strategy is to first collect (if available) and work out the analytical results for Gribov copies 
for the one dimensional lattice and, as far as possible, for the higher dimensional lattices. For 
the one dimensional case, we do have a complete understanding of Gribov copies for the above 
mentioned gauges. For the two dimensional cases, we have been able to derive some lower bounds 
on the number of Gribov copies. Then, we use the recently developed numerical polynomial 
homotopy continuation (NPHC) method to find all Gribov copies for a two dimensional lattice. 
Here, though we could only solve small lattices of the size 3x3 completely, we emphasise that 
this is the first ever result where all Gribov copies are found in more than one dimension. We 
use the NPHC method to study the orbit-dependence of the number of Gribov copies thoroughly 
for these lattices. For the bigger lattices, we move to a traditional method, namely, Monte Carlo 
Conjugate Gradient Method (MCCGM) to attempt to find 'all' minima. 

Our findings are summarised below. 

1. Minimal lattice Landau gauge: The number of minima for the naive functional in one 
dimension with antiperiodic boundary conditions is orbit-independent, but in two dimen- 
sions it is orbit dependent. For periodic boundary conditions, the number of minima is 
orbit-dependent in any dimension. 
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Figure 7: Number of distinct solutions Nsoi found for random orbits on A^ x iV lattices with 
the naive functional. The data points denote the average number of solutions found for each N , 
while the error bars denote the 68% confidence intervals (the 16th and 84th percentiles for the 
number of solutions found). The data for antiperiodic boundary conditions (APBC) are offset 
horizontally for clarity. The N — 8 data are taken from only two random orbits. 

2. Absolute lattice Landau gauge: The number of global minima for any boundary conditions 
for the naive functional is orbit-independent, in both one and two dimensions. This holds 
for both the naive and stereographic gauge fixing functional, in all cases where we have 
obtained results. 

3. Stereographic gauge: We first show that the FP operator for the stereographic gauge is 
generically positive definite, and hence all the stationary points of the corresponding gauge 
fixing functional are minima. We find that the number of minima is orbit-independent in 
one dimension for any boundary conditions, but appears to be orbit-dependent in two di- 
mensions. The number of minima for the stereographic gauge in higher dimensions is much 
higher than even the number of all the stationary points for the naive functional on the 
same lattice. This is exactly opposite to the earlier claim based on the one-dimensional 
results that the number of minima (Gribov copies) is exponentially suppressed for the 
stereographic gauge case compared to the naive gauge. It is worth noting that the stere- 
ographic and naive Landau gauge should become equivalent in the continuum limit and 
that the difference between the two should be due to lattice artefacts. This may still pose 
problems for simulations at realistic couplings, however. 

4. We have not found any clear evidence for the conjecture that the first Gribov region 
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Figure 8: Distribution of the number of distinct minima A^soi found for random orbits with the 
naive functional on different lattice sizes. 

becomes equivalent to the fundamental modular region in the iV — > cxo limit, by studying 
the distribution of functional values at the minima and their median distance from the 
global minimum. It may well be that the convergence to the infinite volume limit is so 
slow that it would not be observed on these small lattices. It is also possible we would 
need to go beyond strong coupling and approach the continuum limit to see this. That is 
beyond the scope of this study. 

Based on our results, we conclude that it would be too cumbersome to fix a gauge using the 
Faddeev-Popov procedure on the lattice with averaging over all Gribov copies, due to the orbit- 
dependence of the Gribov copies in the minimal lattice Landau gauge and stereographic gauge. 
The naive lattice Landau gauge, as is well-known, suffers from the Neuberger 0/0 problem. 
Only the absolute lattice Landau gauge seems to rescue the situation in that it evades the 
Neuberger 0/0 problem by construction and the number of Gribov copies (global minima) is 
orbit independent. However, finding the global minimum of such functions is in general an NP 
hard problem, and it is only if typical local minima become indistinguishable from the global 
minimum that this becomes a viable alternative. Although we find no evidence that the function 
values of the minima as a whole accumulate near the global minimum, we do find that values 
near the global minimum are more likely to be found. It remains an open question, though. 
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Figure 9: Distribution of the values of the naive gauge fixing functional, AF = i^med ~ ^min7 for 
random orbits on different lattice sizes. Left: periodic boundary conditions; right: antiperiodic 
boundary conditions. 
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Figure 10: Number of distinct minima found for the stereographic gauge, for 3 different random 
orbits, as a function of the number of initial guesses iVgucss- Also shown is the result for the 
trivial orbit. 



27 



whether this is still the case for larger lattices and other gauge groups. It is also not known 
whether these near-global minima also exhibit the same qualitative behaviour for quantities such 
as gauge and ghost propagators as the absolute minimum. 

In the context of gauge fixing using a weighted average over Gribov copies as proposed in 
Refs [40 , our results shown in Fig. IH] suggest that the functional values tend to a smooth 
distribution as N grows. This means that a weighted average would be feasible, unlike for 
example the case where there was a gap between the global minimum and all other minima, 
resulting in a potentially serious overlap or ergodicity problem in Monte Carlo simulations. On 
the other hand, if the results in Sec. |4] hold also for larger lattices (and other gauge groups), 
then higher order stationary points are exponentially suppressed compared to the minima, so 
an algorithm that samples only minima (or only some saddle points) may provide a reasonable 
approximation to the full functional integral. 

We should reiterate that Gribov copies in compact U(l) gauge theory are purely a lattice artefact 
with no continuum counterpart, while SU(-/Vc) gauge theories exhibit continuum Gribov copies on 
top of those which are lattice artefacts. However, understanding the nature of Gribov copies in 
compact U(l) is crucial to a full understanding of Gribov copies in SU(A'^c) theories containing 
compact U(l) as a subgroup, and in particular the Neuberger 0/0 problem will be resolved 
once the problem is resolved for compact U(l). Our finding that the stereographic gauge does 
not provide a satisfactory resolution of this problem can therefore be expected to be directly 
applicable to other gauge groups. The bounds on the number of Gribov copies derived in 
Appendix [B] are also valid for all SU(iVc) groups. 
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A Absence of Neuberger problem in stereographic gauge 

Here we sketch the proof of the statement that there is no Neuberger 0/0 problem for the 
steregraphic gauge in any dimension with either APBC or PBC. Essentially, the proof boils 
down to showing that the corresponding FP operator is generically a positive definite matrix 
and hence there is no cancellation among signs of the FP determinants at Gribov copies. 

We first note that the matrix Mpp corresponds to a hessian matrix of a functional with nearest- 
neighbour coupling so it can be decomposed using projection matrices (i.e., matrices with entries 
in { — 1,0,1}) with only two entries per row (or column). Thus, with APBC, the FP operator 
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for the stereographic gauge Eq. ^^ can be decomposed as 



d 



M^pp ^ }^{M,,D^^Mp, (33) 

11=1 

where Df^ are diagonal matrices with the diagonal entries sec^(((/)i^p + 0i+/i — di)/2) with i 
running over all lattice sites and fi running over the lattice dimension d. Here, the matrices M^ 
consist of { — 1,0,1} as their entries, i.e., the i-th row consists of entries corresponding to the 
nearest-neighbour interaction at the i-th site in the /i-th direction. For the APBC case, Af^, are 
non-singular matrices because we have removed the global gauge degree of freedom using the 
boundary condition. Since M ^D'tMj^ are positive definite matrices for all /z = 1, . . . , d, according 
to Sylvester's law of inertia (which states that if A = CBC^ where A is a real symmetric matrix 
and C and B are real matrices with B diagonal, then the number of positive, negative and zero 
eigenvalues of A is the same as those of B), and because the sum of positive definite matrices is 
also a positive definite matrix, Mpp is strictly a generically positive definite matrix. Therefore, 
there is no cancellation of signs of the corresponding FP determinants, and hence, no Neuberger 
0/0 problem. 

With PBC, a decomposition of the corresponding Mpp is more involved so we use another 
method here. We first separate Mpp for different /i in Eq. ([T4|) . i.e., Mpp = J2n{-^Fp)fj. fo'" s- 
d-dimensional square lattice. Now, we consider the quadratic form for the symmetric matrices 
[Mpp)^ with an N'^ x 1 vector y ^ whose elements are y,; where i — {ii, . . . ,id) with ii, . . . ,id 
running over 1, . . . , A^ for the square lattice. Thus, it is straightforward to check that for any /x, 

fiMPp)^. y = ^sec^ -^(2/,+^ - y.)^ (34) 

i 

which is only if all yi are equal, which is the constant zero mode and is strictly positive for 
all other cases. Also, these yi are taken to be following the PBC, in this expression. Thus, the 
matrix (Mpp)^ is positive semi-definite for all /x = 1, . . . ,d. Since the sum of positive semi- 
definite matrices is a positive semi-definite matrix, Mpp is also a positive semi-definite matrix. 
Once we remove the global gauge degree of freedom by taking one of the angles to be zero, say 
(^{N,...,N) = 0, then Mpp is strictly positive definite matrix. Thus there is no Neuberger 0/0 in 
the PBC case either. 



B Lower bounds on the number of Gribov copies 

Here, we show how the Morse theoretical interpretation of the lattice Landau gauge can give 
a lower bound on the number of Gribov copies n. We briefly explain the Morse indices, Betti 
numbers and their relationship with x(M) via the Poincare polynomial P{z) of M along with 
Morse inequalities. We then obtain an expression for the lower bound on n at a given orbit and 
calculate! it expHcitly for compact U(l), SU(2) and SU(A''c). The ultimate goal of this discussion 
is to get a generic (orbit-independent) lower bound on n for the naive and stereographic gauge 
to get a guide on our numerical results, so we ignore the explicit orbit-dcpendcnce here. 



Here, we omit the explicit mention of the orbit-dependence of n. 
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Let Ki be the number of critical points of a height function h{x) with its Hessian at these critical 
points having i negative eigenvalues. Then x(M) = X)i(~l)'^ii where now the sum over i runs 
from to the dimension of M. 

In addition to the Euler characteristic, the Bctti numbers are closely related topological invari- 
ants of a manifold. Firstly, a homology group is a measure of the hole structure of a manifold, 
or more specifically, a topological space. There may be several homology groups of a manifold. 
The ith Betti number hi is defined as the rank of the ith homology group (the reader is referred 
to [55] for details on Betti numbers and their relation to the Euler characteristic) . 

The Poincare polynomial is defined as P{z) = ^^ hiZ^ for an arbitrary real variable z. Then, it 
turns out that 

P(l) - E &. , P(-l) - Y.^-iyh = X(M) , (35) 

i i 

where the sum runs from to the dimension of M. 

There are two types of inequalities, called Morse inequalities, which relate hi and Ki. The weak 
Morse inequality states that 

Y^K,>Y^b, (36) 

i i 

and the strong inequality states that 

K, > h, (37) 

for all i. Due to the weak Morse inequality and using the fact that ^^ Ki is the total number 
of critical points which is n for a lattice Landau gauge fixing functional, we have 

P{l) = Y^b,<Y,K^^n, (38) 

i i 

where for a d-dimensional lattice the sum runs over all the lattice sites, i.e., i = 0, . . . ,N'^ for 
the APBC case and i = 0, . . . , N'^ — 1 for the PBC case. Thus P(l) is a lower bound on n. It 
now remains to calculate P{z) and hence ultimately P{1) for the gauge-group manifold. 

To calculate the corresponding P{z), first note that Px^y{z) = Px{z)Py{z) for a product space 
manifold X xY oi X and Y , where Px{z) and Py{z) are the Poincare polynomials of X and 
Y respectively. For the naive lattice Landau gauge for the compact U(l) case with APBC, the 
corresponding manifold is (S^)^ and the corresponding Poincare polynomial is 



jyO jyO jya 



i 



P^s^).4^) = llPs^{z)^l[{l + z) = {l + zf =^ 

fc=l fc=l i=0 

.P^si)n4^) = 2^'<^i^, = n, (39) 



4=0 



where we have used the fact that Pgi(z) = 1 + z, with the Betti numbers 6o = 1 = ^i a-nd all 
others being zero for S^. Thus, n must be greater than or equal to 2^ in this case. Moreover, 
we can verify that P(— 1) = (1 — 1)^ = = x{{S^)^ ) as expected. Also, we identify the Betti 
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numbers of this iV'-torus as 6^ = ( ^ ), for all i — 0,...,N, by comparing the above equation 
with Eq. ([55)1 . Similarly, for the PBC case, the corresponding manifold is (S^)^ ~^, i.e., 

fc=l fc=l 1=0 ^ ^ 

.•.-P(5i)(«--i,(l) = 2^'-i < ^ i^, - n . (40) 

i=0 

So, the corresponding Betti numbers for {S^)^ ~^ are bi — [ r^) for i — 0, ..., A^'' — 1. 

For SU(2), the group manifold is S^ and the corresponding Poincare polynomial is Ps3{z) = 
(1 + z^). Thus, for the naive lattice Landau gauge with PBC the corresponding Poincare 
polynomial is 

N-l 

^(S3)<«-l)W = 1[{1+Z') = {1 + ZX'-' 

fe=l 
.■.P(,3)(«.-.,(l) - 2^^-i<n^^(2), (41) 

giving the same lower bound for the number of Gribov copies as that of the compact U(l) case. 
For a generic SU(A^c) with Nc > 2, the group manifold is S^ x S^ x ...5^^"="^ and the correspond- 
ing Poincare polynomial is (1 + 2:^)(1 + z^)...(l + z^^~^). Thus, for the naive lattice Landau 
gauge with PBC, 

giving a larger lower bound than that of the compact U(l) and SU(2) cases. 

On the other hand, the corresponding manifold for the stereographic gauge is M.^ for the 
compact U(l) case with APBC. Here, bo — 1, for R, and all other bi = for i = 1, . . . ,N. So 
the corresponding Poincare polynomial, with Pm{z) = 1, is 



pr„.(z) = n(i)=i'^'=i 



fc=i 

.•.P„„.(-l) = Fk«''(1) = 1<"- (43) 

Thus, the lower bound on the number of Gribov copies is exponentially suppressed from 2^ for 
the naive gauge to 1 in the stereographic gauge for compact U(l). Similarly, with PBC, 



^««-.w = n (i)=i^'-'=i 



fe=i 

.■.Pr„._i(-1) = P(R)«.-i(l) = l<n. (44) 

Furthermore, the corresponding Betti numbers of R^^^ are bo = 1 and fe^ = for i = 1, ...,N 
for the APBC case and bo = I and bi = with i = 1, ..., A^^^ - 1 for the PBC case. 
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C Gauge-fixing functional at minima 

Here we want to describe the function values at minima of the gauge-fixing functional in the 
infinite volume limit. 

In one dimension, we have the complete classification of Gribov copies for the naive Landau 
gauge with PBC and odd N, given in Eq. (I^Ul) : the minima satisfy 

(j)1 = (fj^^ = ^ + — mod 27r , ? e {0, . . . , TV - 1} , with ^ - ^*=i ^^ 



together with the condition cos (/)^ > 0. The functional value at the minima is then F^ (0) | minima = 

(1 -COS(/)^). 

By the law of large numbers, since each 0^ is independent and identically distributed, the mean 
value ^ — >■ as A^ ^ oo. Assuming that N is large enough to ignore the 4> term, we notice that 
the stationary points with a random orbit in the thermodynamic limit are the same as those for 
the trivial orbit. Hence in the thermodynamic limit, the first Gribov region of every random 
orbit is identical to the first Gribov region of the trivial orbit. This is a surprising result since 
for finite N, the Gribov regions are orbit dependent. 

In the limit where — ?► 0, the condition cos0^ > becomes equivalent to I < ^ or I > ^. 
Since the function values are the same for minima inZ G {0,...,-^} and I G {^, . . . ,N — 1}, let 
us choose Z e {0, . . . , ^}. Gomputing F — 1 ~ cos{2nl/N) gives us values in the interval [0, 1]. 
We have partitioned [0, 1] into ^ + 1 pieces (the size of the set {0, . . . , ^}). Call the value of 
the I'th piece pi, where pi corresponds to the function value oi I — i. For sufficiently large N, 
the set of all pi will become a countably infinite set obtaining values that appear to map the 
graph of F{x) = 1 — cos(a:) for < a; < -1. 

However, for values oi I '^ N, cos(27r//A^) w 1 and F Ki 0. In our partition of [0, 1] this is 
equivalent to [0,Pk] — >■ as iV ^^ cxd for a particular pk- As N is assumed large, the rest of 
the interval stretches towards zero because we divide by N. Hence the function value of zero is 
hugely degenerate for N sufficiently large with all other function values that are not sufficiently 
near the origin being unique. We now take into account the other interval I ^ {^ , . . . , N — 1} 
which has the same function values as the first interval. Therefore F = is hugely degenerate 
and all other function values not sufficiently close to the origin are 2-fold degenerate. Hence in 
the thermodynamic limit, minima occur taking function values between [0, 1] but zero occurs 
dramatically more often. 

This argument is supported by the numerical results shown in Figure [H We see that even for 
small lattices sizes, the function values at the minima start to accumulate near zero. 
For the naive Landau gauge with PEG in two dimensions, let 6* be a stationary point. Assume 
that this satisfies some conditions Ci,j = 4>i+i.j^x + ^i+i,j ~ ^ij" and di^j = ipij+i^y + dij+i — Oi.j- 
Then the gradient equations become 

sin(ci+ij) - sin(cij) + sm{dij+i) - sin(di,j) = . (45) 

Let us now assume that Cij, dij satisfy the one dimensional gradient equations 

sin((ii,j+i) - sm{dij) = , 

sin(ci+ij) - sin(cij) = . (46) 
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Then these satisfy the 2-diniensional gradient equations, and we have found stationary points 
for which the gradient equations decouple into lower dimensional ones. For a particular orbit, 
we may then find a higher dimensional stationary point by solving the one dimensional gradient 
equations for each hypersurface with j fixed in Cij and i fixed in dij. However, optimising the 
one dimensional gradient equations for Cij with j fixed fixes the j'th column of lattice sites. 
With j G {f, . . . ,N}, this fixes all lattice sites. Therefore, we cannot vary the 9ij in order to 
optimise dij for fixed i with a random orbit. We must impose the trivial orbit. The gradient 
equations are solved only if dij+i — cqu^ij\ where G is some lattice symmetry which imposes 
sin((iij+i) = sm{dij) so that the gradient equations are simultaneously optimised. Now we 
can optimise each hypersurface Cij for fixed j but it is only the subset of these SPs which 
simultaneously imposes the gradient equations of di^ that satisfy the two dimensional gradient 
equations. Using the analytic formula for one dimensional SPs with the trivial orbit from section 
3.1[ we can find these two dimensional SPs. As noted, the minima occur at all qk = and thus 
the two dimensional minima we can construct are 



2ttI _ _ 2ttG{1) 



c^J■ = -irr , dij = cg(^j) = — ^^ , (47) 



where /e{0,...,^,^,...,iV-l}. 

We can explicitly construct some of these minima. With the above conventions, first note that 
requiring dij = {G{1) — 0) simultaneously solves the gradient equations for each I. This is 
easily visualised by first fixing a lattice site to impose boundary conditions and then layering 
a fixed hypersurface of a f-dimensional SP in a constant direction until a 2-dimensional lattice 
is filled. Putting this back into the definition of F it is easy to see that, in this specific case, 
the two dimensional minima has the same function value as the one dimensional minima. The 
Hessian (FP operator) in this G{1) = case is given by 

^FP = c,f L = iM'p%.k6j,i ~ ((5,+i,/ + <5,_i,, - 25,- 05^,fe (48) 

OViAOVkl 



where Af^p is the Hessian of the one dimensional minima, which is a A^ x A^ matrix. From (|48)) 
we see that My^ + 2/ occur in blocks along the main diagonal. The other non zero terms come 
from the second term in (j48p . They are upper and lower diagonal —1 terms starting at the 
(f , A^ + 1) and (A^ + 1, f ) matrix element. Boundary conditions give diagonal — f terms starting 
at the (1, A^^ — A^) and (A^^ — A^, f) matrix element. To fully impose PBC we must remove 
one linearly dependent gradient equation. This reduces the Hessian to a (A^^ — 1) x (A^^ — 1) 
matrix. These minima were observed numerically. More minima can be constructed by defining 
G and we can find others by the permutation symmetry as described in section [3.21 Therefore 
we can construct multiple sets of two dimensional minima were each set is in a one-to-one 
correspondence with the set of all one dimensional minima. 

Since we can construct a number of 2-dimensional minima from each one dimensional minimum, 
the two dimensional first Gribov region will contain multiple analogues of the one dimensional 
first Gribov region. By the above arguments, since most function values of the 1-dimensional 
minima accumulate at zero in the thermodynamic limit, this behaviour must be imitated locally 
by the 2-dimensional first Gribov region. As this construction of minima can be extended to 
higher dimensions, this argument generalises. 
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We cannot comment on the global behaviour of the two dimensional first Gribov region as an 
analytic classification of all two dimensional minima remains elusive. However, the numerical 
results from Figure |B] show that the function values of 2-dimensional minima with a trivial 
orbit appear to decay towards zero, similar to the one dimensional case. This suggests that the 
imitation may be global for the trivial orbit. Since this construction only works for the trivial 
orbit, there is no reason to suggest that the two dimensional first Gribov region for a random 
orbit behaves locally like the one dimensional first Gribov region. Indeed, this difference between 
the trivial and random orbit case was observed in section 15.2.11 and particularly in Table 21 
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N 


3 4 5 6 


7 


Norh 


100 100 100 50 


10 


-^*gucss 


2 • 10"* 5-104 2 • 10^ 2 • 10^ 


6-107 


Periodic boundary conditions 


(Nsoi) 


1.9 4.5 14 70 


373 


68% CI 


1-3 2-7 7-21 40-111 


276-534 


\-^inin/ 


OMtll 0.48±i 0.44±i 0.44±? 


0.43±l 


\-f^mcd -t^nun) 


0.12-12 0.11-8 0.15-7 0.14-4 


0.15-3 


Antipcriodic boundary conditions 


{Nsoi/2) 


2.7 7 20 93 


717 


68% CI 


2-4 4-10 12-29 42-165 


284-1493 


\-^min/ 


0.43-11 0.43-8 0.41-7 0.40-4 


0.42-4 


\-^mcd -^min/ 


+21 +9 +5 +4 

0.17-11 0.15-9 0.14-4 0.14-3 


0.13±l 



Table 4: Number of random orbits A^orb and initial guesses iVguoss used for the different lattice 
sizes using the naive functional, together with results for periodic and antipcriodic boundary 
conditions. Agoi denotes the number of minima found and 68% CI is the 68% confidence interval 
in the number of minima. -Fmin is the absolute minimum of the gauge fixing functional, and 
its median value at the local minima is -Fmod- The errors in the last digit (s) denote the 68% 
confidence interval. 
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